import numpy as np

def update_linear_velocity(current_acceleration, current_velocity, time_interval):
    """
    Update linear velocity given the current acceleration and current velocity.
    
    Parameters:
    current_acceleration: array-like, shape (3,)
        An array containing the current acceleration in m/s^2 for each axis [x, y, z].
    current_velocity: array-like, shape (3,)
        An array containing the current velocity in m/s for each axis [vx, vy, vz].
    sampling_frequency: int
        Sampling frequency of the accelerometer in Hz.
        
    Returns:
    updated_velocity: array-like, shape (3,)
        An array containing the updated velocity in m/s for each axis [vx, vy, vz].
    """
    

    
    # Update velocity using Euler method
    updated_velocity = current_velocity + current_acceleration * time_interval
    
    return updated_velocity

# Example usage:
# Let's assume we have an initial velocity
current_velocity = np.array([0.0, 0.0, 0.0])
sampling_frequency=60
time_interval = 1 / sampling_frequency
# Let's simulate receiving new acceleration data over time
for _ in range(100):  # Assume 100 new data points
    # Simulate new acceleration data
    current_acceleration = np.random.normal(loc=0.0, scale=0.1, size=3)
    
    # Update velocity
    current_velocity = update_linear_velocity(current_acceleration, current_velocity, time_interval)
    
    print(f"Updated Velocity: {current_velocity}")